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Abstract. The expected metal enrichment of the intra-cluster medium (ICM) 
and the partition of metals between cluster galaxies and the hot ICM depends 
on the stellar Initial Mass Function (IMF). The choice of the IMF in simulations 
of clusters has also important consequences on the "cold fraction", which is a 
fundamental constraint on cluster physics. 

We discuss the chemical enrichment and the cold fraction in clusters as predicted 
with different IMFs, by means of a straightforward approach that is largely in- 
dependent of the details of chemical evolution models or simulations. We suggest 
this simple approach as a guideline to select the input parameters and interpret 
the results of more complex models and hydrodynamical simulations. 
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1. Introduction 

The chemical enrichment of the intra- 
cluster medium (ICM) has been extensively 
discussed in literature by means of chemi- 
cal evolution models of elliptical galaxies 
with galactic winds, and a variety of sce- 
narios has been advanced (see the review 
by Matteucci, these proceedings). Very re- 
cently, cosmological hydro-dynamical sim- 
ulations of cluster formation have been de- 
veloped, that can follow self-consistently 
the star formation and chemical enrichment 
history of cluster galaxies and of the ICM 
(Valdarnini 2003; Tornatore et al. 2004; 
Romeo et al. 2004). Star formation has 
important effects on the hydro-dynamical 
evolution of the cluster, via energy feedback 
from supernova? and metal enrichment of 
the ICM: the first effect contrasts, the sec- 



ond one boosts the cool-out of the hot gas. 
The chemical evolution of the ICM is not 
only an interesting issue per se, but an im- 
portant ingredient of the global physical 
evolution of clusters. As a consequence, the 
choice of input parameters and "recipes" 
related to star formation (notably, the stel- 
lar Initial Mass Function and the imple- 
mentation of sub-grid feedback effects) is 
crucial for the results of the simulations. 

In this paper we outline a simple proce- 
dure to estimate the chemical enrichment 
of the ICM, the partition of the metals be- 
tween stars and ICM, and the cold frac- 
tion expected after an assumed Initial Mass 
Function (IMF). This provides a guide- 
line to select the optimal input parame- 
ters of the simulation; and to distinguish, 
in the results of a complex and fully self- 
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Fig. 2. Mass fraction R of returned 
gas from dying stars and complementary 
locked-up fraction M* for a Salpeter SSP. 



Fig. 1. (a) Rates of SN II and SN la (num- 
ber per year) for a 1 M Q Salpeter SSP. (b) 
Cumulative iron production from SN II and 
SN la. (c) Luminosity evolution of the SSP. 
(d) Evolution of IMLRssp. 

consistent simulation, what is merely a con- 
sequence of the adopted IMF, and what is 
an effect of the interplay between hydrody- 
namical evolution and star formation. 

2. IMF, metal production and 
partition between stars and ICM 

Let's choose an IMF and consider a burst of 
star formation (a Single Stellar Population, 
SSP) with stellar masses distributed ac- 
cordingly; we can compute the expected 
rates of SN II and SN la, the correspond- 
ing rate of production of metals (for ex- 
ample, iron: Mp e (t)), and the luminosity 
evolution of the SSP Ls{t). We can then 
define the Iron Mass-to-Light ratio typical 
of that IMF as IMLR SSP (i) = Fig.ffl 
illustrates this procedure for the Salpeter 
(1955) IMF with mass limits [0.1-100] M Q ; 
the detailed calculations can be found in 
Portinari et al. (2004, hereafter PMCS). 

Part of the iron produced will be "eaten 
up" by subsequent star formation episodes, 
to build up the stellar metallicities observed 
in cluster galaxies; only the remaining frac- 
tion is available to enrich the ICM. We need 



to estimate how the total iron produced 
gets partitioned between the stellar popu- 
lations and the ICM (Renzini et al. 1993). 

We can estimate the amount of iron 
locked in the stars as Mp e .* = Zp e .* x M*. 
As to the stellar metallicity Zp e ,*, the star 
mass in clusters is dominated by massive el- 
lipticals, with global metallicities between 
(solar) and +0.2 dex, and [a/Fe] ra- 
tios around +0.2 dex (PMCS and refer- 
ences therein) . The mass in stars M* is not 
directly observable; it is usually inferred 
from the global luminosity of cluster galax- 
ies, but the conversion factor, the stellar 
Mass-to-Light ratio (M*/L), depends on 
the IMF. Thus, the mass M* locked in stars 
is best computed self-consistently from the 
assumed IMF (Fig.|3Jl; M* drives the actual 
amount of iron Mp e ^ required to reproduce 
the observed metallicity Zp e ,*- 

We are now able to split the character- 
istic IMLRssp of the assumed IMF into 
the fraction IMLR* = M f'-* locked in the 
stellar populations, and the remaining frac- 
tion IMLR/cm that can be compared to 
ICM observations (Fig. red line vs. red 
shaded area). Notice that IMLR/cm as 
computed here is an upper limit to the ex- 
pected enrichment of the ICM, since all the 
iron not locked in the stars is assumed to 
be expelled from the galaxies. Whether this 
actually occurs, depends on the efficiency of 
extraction/ejection mechanisms (feedback, 
ram pressure, etc.). 
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Fig. 3. (a) IMLRssp for the Salpeter IMF 
split into its components: IMLR* locked in 
the stars and IMLR/cm in the ICM; the 
red shaded area is the observed IMLR in 
the ICM (Finoguenov et al. 2000, 2001). 
(b) Same as (a), for the Silicon Mass-to- 
Light ratio, (c) [Si/Fc] ratio in the stars 
(assumed a priori) and in the ICM. 

The same exercise, described above for 
iron, can be repeated for any other ele- 
ment; in particular for the a-elements that 
directly trace SN II enrichment and star 
formation. Among these, silicon is the best 
measured one in the ICM (Fig. 03). 

From Fig. |3 our simple argument shows 
that with the Salpeter IMF one expects: 

1. equipartition of iron between stars and 
ICM (Fig.EK; Renzini et al. 1993); 

2. an IMLR in the ICM compatible with 
observations, provided the bulk of the 
stellar population in clusters is old (say, 
older than 10 Gyr; Fig. Et); 

3. an uneven distribution (i.e. no equipar- 
tition) of a-elements, that are mostly 
contained in the stars (Fig. Et>) ; 

4. a "chemical asymmetry" between stars 
and ICM (Fig. Et; Matteucci & 



Vettolani 1988; Renzini et al. 1993; 
Pipino et al. 2002): the stars have su- 
persolar [a/Fe] abundance ratios — as 
we assumed a priori after observational 
results — while the ICM has undersolar 
[a/Fe]; the latter result is at odds with 
observations (Finoguenov et al. 2000, 
2001; Baumgartner et al. 2004); 
5. a SiMLR lower than observed (Fig.Efc). 

In a hydro-dynamical simulation adopting 
the Salpeter IMF, these basic predictions 
can be used to check and interpret the re- 
sults. For instance, if more iron is contained 
in the stars than in the ICM, it means 
that the simulated stellar iron abundances 
are in excess of the observed ones and a 
more efficient ejection of iron into the ICM 
needs to be implemented. If the IMLR in 
the simulated ICM is lower than observed, 
either not enough iron has been ejected 
into the ICM, or the stellar populations 
are too young (high luminosity); while a 
low SiMLR and an undersolar [Si/Fe] ratio 
in the ICM are expected with the Salpeter 
IMF and are not a fault of the simulations. 

The results listed above refer to the 
Salpeter IMF, but the procedure can be ap- 
plied to any other IMF; the expected metal 
production and partition between stars and 
ICM will be different. Fig.EJshows the case 
for a typical Solar Neighbourhood IMF, the 
Kroupa (1998) IMF. With a steep Scalo 
slope at the high-mass end, the Kroupa 
IMF has a lower iron production than 
Salpeter (Fig. 0t). Being also "bottom- 
light" with respect to Salpeter (namely, 
with less mass locked in low-mass stars) 
it has a higher luminosity as well (Fig.QJ)). 
All in all, its characteristic IMLR55P is a 
factor of 2 lower than that of Salpeter, and 
it cannot possibly match the observed iron 
enrichment in the ICM (Fig.Et; see PMCS 
for details). Even worse is the outcome 
for the a-elements (Fig. QJl). Henceforth, 
one can predict beforehand that adopting 
a "standard" Solar Neighbourhood IMF 
in cluster simulations can never lead to a 
satisfactory ICM enrichment. Notice also 
that with the Kroupa IMF there is no 
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Fig. 4. (a) Iron production from a Kroupa 
SSP, compared to a Salpeter SSP. (b) 
Luminosity evolution of a Kroupa vs. 
Salpeter SSP. (c) IMLR SS p for Kroupa 
and Salpeter; the Kroupa IMLR55P is split 
into its stellar and ICM components, and 
the latter is compared to observations (red 
shaded area), (d) Same as (c), for the 
SiMLR. (d) [Si/Fe] ratio in the stars and 
in the ICM, for the Kroupa IMF. 



equipartition of iron between stars and 
ICM (Fig. UJ;); iron equipartition holds 
specifically for the Salpeter IMF. 

Fig. |5] shows the expected production 
and partition of elements between stars and 
ICM, for the top-heavy Arimoto & Yoshii 
(1987) IMF (power-law exponent -1.0 vs. 
-1.35 for Salpeter; same mass limits [0.1- 
100] M Q ). For this IMF, more metals are 
available to enrich the ICM than locked 
in the stars, and the IMLR and SiMLR 
in the ICM can match the observed lev- 
els for stellar populations of 5-10 Gyrs of 
age. The expected [a/Fe] ratios in the ICM 
are supersolar, in agreement with observa- 
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Fig. 5. IMLR, SiMLR and [Si/Fe] ratio in 
the stars and in the ICM, expected from 
the Arimoto & Yoshii IMF. 



tions. Therefore, a top-heavy IMF like this 
can potentially lead to satisfactory results; 
whether this is achieved in an actual sim- 
ulation, depends on the simulated star for- 
mation history and ejection mechanisms of 
metals from the galaxies. 

3. IMF, M*/L ratio and cold fraction 

In the previous section we showed that 
a self-consistent account of the mass and 
metals locked in the stars, related to the 
M*/L corresponding to the adopted IMF, 
is fundamental to predict the chemical en- 
richment of the ICM. The M*/Z has also 
important consequences on the cold frac- 
tion, defined as the ratio M*/Micm be- 
tween stellar mass and hot ICM gas mass 
(the cold gas component has a negligi- 
ble mass contribution). The cold fraction 
is crucial to distinguish, for instance, the 
mechanisms responsible for the "entropy 
floor" in low temperature clusters. The 
larger the cold fraction, the larger the role 
played by galaxy formation that removes 
low entropy gas; with a small cold fraction, 
instead, strong supernova (pre)heating or 
other feedback sources are needed to set 
the entropy floor (Bryan 2000; Balogh et al. 
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2001; Valdarnini 2003; Tornatore et al. 
2003). Consequently, a correct estimate of 
the cold fraction is a major constraint for 
scenarios of cluster evolution. 

While Micm can be inferred directly 
from the observed X-ray emission, the stel- 
lar mass can be estimated only indirectly 
from the observed luminosity, typically in 
the B-band, via an assumed M*/L. In liter- 
ature, very discrepant values can be found 
for the cold fraction — or its reciprocal: 

Micm 



a 

O0.2 - 



O0.1 
o 



= 2-20 



This discrepancy is mostly due to the as- 
sumed M^/L; the determinations of the di- 



rectly observed quantity 



ICM 



(listed be- 



low as compiled by Moretti et al. 2003) 
are in much better agreement, within 30% 
around: 



M. 



ICM 



L, 



= 30 h 



.1M0 



36- 



(h = 0.7) 



Micm 
Lb 



h 2 —Si 



30 David et al. (1990), Hydra A 

30 (± 13) Arnaud et al. (1992) 

30 White et al. (1993), Coma 

19 Cirimele et al. (1997) 

31-44 Roussel et al. (2000) 

>21 Finoguenov et al. (2003) 



With this typical value of 



M 



ICM 



, the cold 



fraction depends on the assumed M„/Lb 
as shown in Fig. El Displayed are the values 
of M„/L B assumed by White et al. (1993) 
and Balogh et al. (2001), by Fukugita et 
al. (1998; averaged over the typical mor- 
phological mixture of cluster galaxies out 
to the Abell radius, or for spheroids only); 
and the M*/Lb values corresponding to 
the Salpeter, Kroupa and Arimoto & Yoshii 
IMF, assuming an age of 10 Gyr for the 
stellar populations. The estimated cold 
fraction ranges between 10 and 20%; while 
such factor of 2 is of little relevance for the 
"baryon budget" addressed by White et al. 
and Fukugita et al., since the baryonic mass 
in clusters is anyways largely dominated 
by the hot ICM, it has a drastic impact 
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Fig. 6. Cold fraction M^/Micm in clusters 
of galaxies as a function of the assumed 
M^/Lb ratio, for a typical observed value 
of Micm /L B =36 M /i . 

on chemical enrichment, cold fraction and 
cluster modelling. Fig. shows that, if one 
assumes a Salpeter (or Arimoto & Yoshii) 
IMF, the observed M ^ M corresponds to a 
cold fraction around 20% for suitable ages 
of the stars in cluster galaxies (~10 Gyr, 
suited to reproduce the IMLR: Fig. 0t)- 
The widely quoted cold fraction of 10% can 
be reached only by simulations that assume 
an IMF corresponding to M^/Lb values as 
low as assumed by e.g. White et al. (1993). 

The correct constraint for a simulation 
is the cold fraction corresponding to the as- 

MiCM 



sumed IMF, based on the observed 
Equivalently, for the simulated cluster one 
can compute self-consistently the luminos- 
ity relevant to the assumed IMF and com- 
pare directly to the observed M [ C B M ■ 

Alternatively, if simulations aim at 
reaching a cold fraction of 10%, one should 
select beforehand an IMF (quite "bot- 
tom light") able to reach as low values of 
M*/Lb as assumed by White et al., for the 
old ages typical of elliptical galaxies. 

4. Conclusions 

We remark the following crucial points for 
modelling the cosmological evolution and 
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chemical enrichment of clusters of galaxies 
self-consistently. 

— The amount of mass and metals locked 
in the stellar component is not neces- 
sarily negligible, depending on the as- 
sumed IMF and corresponding M*/L. 

— The observed metal content in stars can 
be used as a constraint for the effi- 
ciency of feed-back and of metal dis- 
persion in the simulations: the wind 
efficiency cannot be indefinitivcly en- 
hanced, because the amount of metals 
in the stars should also be accounted 
for, and a sizable fraction of the metals 
produced may thus not be available to 
enrich the ICM. (Though the problem 
of present-day simulations is rather the 
opposite, namely to lock too much met- 
als in the stars; Tornatore et al. 2004; 
Romeo et al. 2004). 

— Once an IMF is chosen for the simu- 
lations (preferably a bottom-light IMF 
with high-mass slope shallower than 
Scalo, see PMCS), we suggest to com- 
pute the corresponding partition of 
metals between stars and ICM with the 
simple procedure outlined in this pa- 
per. Such expected partition can then 
be used to test the numerical results. 

— The crucial constraint of the cold frac- 
tion should be consistent with the 
M*/L relevant to the assumed IMF. As 
the observational quantity is luminosity, 
we suggest to compute the luminosity of 
the simulated cluster consistently with 
the adopted IMF, and compare it di- 
rectly to the observed Micm/L. Red or 
IR luminosities are preferred, since they 
probe the stellar mass with a lower sen- 
sitivity to the age of the stellar popula- 
tions and to recent star formation activ- 
ity; Micm /Lk estimates are becoming 
available (Lin et al. 2003). 
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